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The Thermal Infrared Sensor (TIRS) 1 


• Purpose: 

- One of two instruments on the Landsat Data Continuity 
Mission (LDCM) 

- Provide 100-m ground-resolution thermal imagery 

- Ensure continuity with previous Landsat thermal instruments 

• Design: 

- Push-broom type imager 

- Refractive, 4-element telescope 

- Operates in two thermal bands: 10.9 and 12.0 /vm 

- Three Quantum Well Infrared Photodetector 
(QWIP) 2 arrays 

1 D. Reuter, et al., “The Thermal Infrared Sensor on the Landsat Data Continuity Mission,” in Geoscience and Remote Sensing Symposium (IGARSS), 
2010 IEEE International, pp. 754-757 (2010). 

2 M. Jhabvala, et al., “QWIP-based thermal infrared sensor for the Landsat Data Continuity Mission,” Infrared Physics & Technology 52, 424-429 (2009). 



TIRS Focal Plane 



• Consists of three QWIP arrays 

- Each QWIP array contains 2 spectral bands, centered on 
10.9 and 12.0 /vm 
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Purpose of Test 



• Each QWIP array can be displaced from the nominal focal 
plane in 3 degrees-of-freedom: 

- Tip about the cross-track axis 

- Tilt about the along-track axis 

- Translation along the optical axis 

• No active focus mechanism in the instrument 

- Focal plane must be “Set & Forget” at system integration 

- Some small amount of re-focusing allowed on orbit by controlling 
telescope temperature 

• Need an accurate method of determining system focus for a 
given focal plane placement 

- Manufacture & install shims to correct 
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Methods 



• Point-spread function (PSF) focal-sweep data was collected at 18 
points across the field-of-view (3 in each band on each QWIP) 

- Defocus achieved by sweeping the object point through a limited range along 
the optical axis 
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• Each focal sweep was analyzed to determine the nominal defocus of 
the detector array at that field point 

- Defocus at each field point is converted to tip, tilt and displacement of each 
QWIP array 

• Two methods for determining defocus: 

- Image-based wavefront sensing (phase retrieval) 

- Gaussian-fit model prediction 


Methods: Phase Retrieval 

• Employs a nonlinear optimization-based algorithm 3 

• Challenging phase-retrieval problem: 

- Undersampled data, Q = A(F/#)/( pixel size) = 0.78 at A = 12.0 /vm 

- Bandwidth (AA / A) * 7% 

- Limited diversity defocus of »±0 .7 A at 12.0 jum 

- Barely unresolved target: object pinhole diameter = 1 detector pixel 

• Jointly optimizes the pupil phase using all the images in the 
focal sweep, allowing diversity to vary for each image 

- Diversity defocus is implemented as varying amounts of 
Zernike defocus (Z 5 ) 

3 S.T. Thurman and J.R. Fienup, “Complex pupil retrieval with undersampled data,” J. Opt. Soc. Am. A 26 , 2640-2647 (2009). 



Methods: Phase Retrieval 




Object Position [mm] 


Methods: Gaussian Fit 
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Methods: Gaussian Fit 



• Plot the Gaussian a fit for each PSF against source 
defocus position: 
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Data: Single Field Point 

12.0 /jm Data 
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Results: Estimated PSFs 

12.0 jiym Data 



NRMSE = 0.1183 


NRMSE = 0.1797 


NRMSE = 0.1643 
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Results: Comparison 

10.9 jiym Data 
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RMS Error = 0.0105 mm 
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Conclusion 



• Phase retrieval was used to estimate defocus at various 
points across the field-of-view 


• Challenging wavefront-sensing problem due to undersampled 
data, minimal diversity, and nearly-resolved object 


• Phase-retrieval results were corroborated by an independent 
Gaussian-fit technique 


• Focal-plane displacements from the two techniques agreed to 
within an RMS error of 

- 0.0105 mm for A = 10.9 ym data 

- 0.0101 mm for A = 12.0 /jm data 



BACKUP 
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Definition of “Best Focus” 

in the Presence of Aberration 



Minimum WFE 
(Zernike Defocus = 0) 



Minimum PV Spot Size Minimum RMS Spot Size 

(Circle of Least Confusion) (Gaussian Min.) 
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1. Introduction 

Image-based wavefront sensing techniques, such as phase-diverse phase retrieval (PDPR), are becoming 
increasingly common tools for aligning complicated optical systems, as well as for performing optical 
metrology at the component level [1, 2]. Generally speaking, image-based techniques are more 
straightforward to apply when the detected point-spread function (PSF) intensity data is critically sampled 
(Q = 2, where Q = 2(F/#) / (pixel size)). Recent work has further improved the performance of image- 
based techniques when applied to undersampled data [3,4], specifically when Q < 1. 

Building from this earlier work, we have recently performed phase retrieval on under-sampled data 
collected from NASA’s TIRS optical system [5]. We applied the approach of Ref. [3] to estimate the low- 
order wavefront error of TIRS across the FOV. Of particular interest is the defocus term of the wavefront 
as it correlates directly to the relative displacement of the detector along the optical axis at each field point. 

2. System Description 

The system under test (SUT) consisted of a 4-element refractive optical telescope and three Quantum Well 
Infrared Photodetector (QWIP) arrays, operating in a 1-pm bandwidth, centered at 12.0 pm (A2/2 = 0.08) 
with a sampling Q value of 0.78. Blackbody point-source illumination was collimated by an off-axis 
parabola (OAP) and injected into the SUT. Scanning the point source across the FOV was accomplished 
by a two-axis steering mirror after the OAP collimator, while introducing defocus was accomplished by 
moving the point source object along the optical axis about the focal point of the OAP. 

In total, 9 field points were examined, the locations of which were determined by active regions on the 
infrared detector [7]. The range of motion available to the point-source object limited the total amount of 
defocus that could be introduced to approximately ± 0.7 2 of Zernike defocus at 12.0 pm. Thus, with the 
challenges of undersampled PSFs and small defocus diversity, data collected from this imaging system 
poses a difficult estimation problem. 

3. Results 

A nonlinear optimization phase-retrieval algorithm processed data from all 9 field-points, using 11 
diversity defocus images at each field point. The global defocus term was recovered by plotting the 
Zernike defocus coefficient for each of the diversity images against the corresponding position of the point 
source relative to the OAP focal point. The zero crossing of a linear fit thus indicates the position at which 
the point source should be focused onto the detector, and therefore how far the detector is from the focal 
plane of the optical system. 

The detector displacement values obtained using the phase-retrieval algorithm are compared to values 
obtained from a second, independent method, which was developed to determine shim adjustments for the 
focal plane assembly. In the alternate method, the detected PSFs were fit to 2-dimensional Gaussian 
functions, and the o of the Gaussian fit was plotted against the point-source object’s position relative to the 
OAP focal point. A parabolic fit to this data then determined the source position that gave the minimum o 
value. This position for each field point was transferred to an as-built optical model of the SUT and OAP 
collimator, which allowed for an estimation of the tip, tilt, and defocus of the focal plane assembly. Setting 



the source defocus back to 0 in this optical model (i.e. collimated input light) produced an estimate of the 
defocus at the detector that could be directly compared to the PDPR estimate. 

Figure 1 shows a plot of the estimated defocus value for each field point, using both the PDPR 
technique and the Gaussian fit method. The two techniques are in good agreement across the FOV, with a 
root-mean- squared difference of 0.010 mm. 
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Figure 1 - Comparison of PDPR results with Gaussian fit technique at each of the nine field points. 



4. Summary 

We present results in which PDPR was used on an undersampled, thermal infrared imaging system with 
limited diversity defocus to estimate the low-order wavefront error terms across the FOV. The estimated 
global defocus was used to determine focal-plane displacement and agrees with an independent Gaussian 
fit technique to within 0.010 mm RMS. 
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